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Abstract 

In  transportation  systems  today,  there  is  a  need  to 
predict  where  a  vehicle  will  be  at  a  given  time  in  or¬ 
der  to  ensure  safety,  expediency  and  efficiency  of  traf¬ 
fic  movement.  There  is  generally  a  plan  of  travel,  but 
outside  forces  (e.g.,  wind  forecasting  error,  navigation 
system  error)  cause  the  actual  path  that  is  followed  to 
be  somewhat  different  from  the  planned  path.  The  path 
of  a  vehicle  is  represented  as  a  vector-valued  curve  in 
three-space.  The  construction  of  the  confidence  region 
about  the  curve  takes  advantage  of  an  assumption  that 
the  deviation  of  the  actual  path  from  the  predicted  path 
will  satisfy  the  conditions  for  a  conditioned  Brownian 
motion  process.  Using  a  cubic  spline  to  estimate  the 
predicted  path,  it  is  possible  to  obtain  parameter  values 
for  the  conditioned  Brownian  motion  process,  as  well  as 
error  bounds  for  constructing  the  confidence  region.  An 
example  is  given  to  illustrate  the  forecasting  technique, 
showing  good  results  in  predicting  the  path,  construct¬ 
ing  the  confidence  region,  and  detecting  when  the  actual 
distribution  of  the  deviation  differs  from  the  estimated 
distribution. 

1.  Introduction 

There  is  a  need  to  predict  where  a  vehicle  will  be  at 
a  given  time  in  order  to  ensure  safety,  expediency  and 
efficiency  of  traffic  movement.  The  location  of  the  vehi¬ 
cle  can  be  described  in  four  dimensions:  lateral  (speci¬ 
fied  in  latitude  and  longitude),  vertical  (a  distance  above 
ground  level  or  sea  level),  and  longitudinal  (time).  The 
lateral  dimension  is  here  expressed  as  a  composite  of 
two  dimensions.  In  rail  systems,  the  lateral  and  vertical 
dimensions  are  easily  predicted,  leaving  only  time  un¬ 
certainty.  In  automotive  systems  (e.g.,  busses),  vertical 
dimensions  are  uninteresting,  lateral  dimensions  c^ln  be 
well-defined  (bus  routes)  or  not  (personal  automobiles), 
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and  time  is  uncertain.  In  certain  amphibian  systems 
(e.g.,  ferries),  there  is  more  uncertainty  in  the  lateral  and 
time  dimensions,  but  still  the  vertical  dimension  is  un¬ 
interesting.  However,  in  underwater  or  air  travel,  there 
can  be  large  uncertainty  in  all  four  dimensions.  In  these 
cases,  there  is  generally  a  plan  of  travel  which  is  fol¬ 
lowed,  but  outside  forces  (wind  forecasting  error,  naviga¬ 
tion  system  error)  cause  the  actual  path  that  is  followed 
to  be  somewhat  different  from  the  planned  path. 

Path  prediction  has  tended  to  emphasize  short-term, 
or  tactical,  path  prediction  of  one  vehicle  at  a  time,  usu¬ 
ally  in  the  arena  of  military  intelligence  (e.g.,  the  predic¬ 
tion  of  the  near-term  positions  of  a  submarine  belonging 
to  the  opposing  force).  The  position  of  an  aircraft  in  the 
air  traffic  control  system  is  predicted  for  two  to  three 
minutes  in  the  future  and  checked  against  the  short-term 
predicted  future  positions  of  other  aircraft  in  the  vicin¬ 
ity  in  order  to  determine  if  separation  minima  are  likely 
to  be  violated.  The  short-term  predicted  position  of  the 
aircraft  is  based  upon  its  reported  current  position  and 
on  a  small  number  of  immediately-previous  reports,  plus 
knowledge  of  whether  the  aircraft  is  near  a  planned  route 
bend  and  the  direction  of  turn  (left  or  right).  The  long¬ 
term  prediction  is  made  only  for  a  limited  number  of  se¬ 
lected  positions  on  the  planned  route  and  is  made  based 
on  the  planned  route  and  the  planned  cruising  speed. 
In  tactical  path  prediction,  past  reports  can  be  used  to 
improve  the  prediction,  by  the  use  of  techniques  such  as 
a  Kalman  filter.  The  extended  Kalman  filter  and  the 
adaptive  Kalman  filter  have  also  been  used  as  a  means 
of  updating  the  predicted  path  of  a  vehicle,  especially  in 
the  design  of  control  systems.  (See  Lefas  and  Thomas, 
1981,  and  Brown,  1983.)  These  filters  are  especially  ap¬ 
propriate  when  a  vehicle  has  an  extended  mission  over 
several  months  and  when  initial  uncertainties  and  sparse 
measurements  are  not  critical  factors  (Brown,  1983,  pp. 
300-307). 

In  this  paper,  the  path  of  a  vehicle  is  represented  as 
a  vector-valued  function.  A  confidence  region  is  con¬ 
structed  about  the  planned  path,  based  upon  previous 


traversals  of  the  same  path  by  the  same  vehicle  or  other 
similar  vehicles,  or  upon  assumptions  about  error  limits. 
The  construction  of  the  confidence  region  takes  advan¬ 
tage  of  an  assumption  that  the  deviation  of  the  actual 
path  from  the  planned  path  will  satisfy  the  conditions 
for  a  conditioned  Brownian  motion  process.  Using  a  cu¬ 
bic  spline  to  estimate  the  predicted  path,  it  is  possible 
to  obtain  parameter  values  for  the  conditioned  Brownian 
motion  process,  as  well  as  error  bounds  for  constructing 
the  confidence  region.  An  example  is  given  to  illustrate 
the  forecasting  technique,  showing  good  results  in  pre¬ 
dicting  the  path,  constructing  the  confidence  region,  and 
detecting  when  the  actual  distribution  of  the  deviation 
differs  from  the  estimated  distribution. 


2.  Path  Modeling 


We  introduce  the  following  terminology  for  the  move¬ 
ment  of  the  object.  The  planned  path  of  the  object  rep¬ 
resents  the  idealized  planned  movement  of  the  object 
in  the  future,  defined  in  terms  of  points  of  trmsition 
in  each  spatial  dimension.  An  actual  path  of  the  object 
represents  the  actual  positions  of  the  object  as  measured 
during  one  realization  of  the  object’s  transversal  of  the 
planned  path.  The  mean  path  of  the  object  represents 
the  movement  of  the  object  “averaged”  over  many  at¬ 
tempts  of  the  object  to  follow  the  planned  path,  i.e.,  a 
mean  of  actual  paths.  The  predicted  path  of  the  object 
is  the  path  that  the  object  might  be  expected  to  follow 
during  any  actueJ  realization  of  the  planned  path. 

We  wish  to  model  the  predicted  path.  In  doing  so,  we 
view  the  error  in  prediction  as  a  Brownian  motion  pro¬ 
cess  with  known  points  of  constriction  corresponding  to 
locations  which  must  be  traversed  (e.g.,  turns).  At  the 
points  of  constriction,  the  error  is  expected  to  be  very 
small.  Thus,  we  model  the  error  between  constriction 
points  as  a  Brownian  bridge.  We  show  that  the  pre¬ 
dicted  path  can  be  estimated  using  a  vector  spline.  It 
is  desired  to  determine  a  confidence  region  or  envelope 
about  the  predicted  path.  The  confidence  region  repre¬ 
sents  a  volume  determined  as  a  function  of  the  error  in 
prediction  and  within  which  the  object  should  be  located 
with  a  given  degree  of  probability. 

We  represent  the  planned  path  as  a  vector-valued 
function  P{t)  =  (pi(i), . . .  ,p„(t))^  of  the  single  real 
variable  t,  for  0  <  ^  <  T.  We  are  typically  interested 
in  the  cases  when  the  dimension  n  is  either  two  or  three, 
and  when  the  planned  path  represents  an  attempt  to  fol¬ 
low  a  straight  line  between  adjacent  pairs  of  points  (as  in 
the  movement  of  an  aircraft).  In  these  cases,  when  the 
value  of  P  is  known  for  m  distinct  values  ti, . . . ,  >  the 


planned  path  can  be  visualized  as  a  connected  sequence 
of  straight  line  segments  -  the  vector- valued  function 


P(0  = 


t-tj 
1  —  fi 


^•+1  —  ^ 
^•+1  ft 


(Pm 


for  t  €  (f,-,f,+i]  and  0  <  ti  <  . . .  <  tm  <  T  connected 
at  the  points  P{ti)  (called  joints).  It  can  be  assumed 
(although  this  is  not  essential)  that  each  joint  represents 
a  change  in  direction  in  one  of  the  spatial  dimensions. 
Because  of  the  abrupt  changes  in  direction  at  the  joints, 
the  planned  path  belongs  at  best  to  C®. 

An  actual  path  Y  (f )  represents  the  locations  traversed 
by  the  object  while  trying  to  follow  the  planned  path. 
The  locations  observed  on  the  actual  path  are  sampled 
values. 

The  mean  path,  also  a  vector-valued  function,  rep¬ 
resents  the  mean  of  a  large  number  of  realizations  of 
the  actual  path.  The  mean  path  m{t)  is  m{t)  =  E[Yt], 
Each  meein  path  is  assumed  to  be  smooth  in  the  sense 
that  any  functional  representation  should  have  continu¬ 
ous  first  and  second  derivatives.  That  is,  we  assume  any 
vehicle  will  move  in  such  a  manner  that  velocity  and  ac¬ 
celeration  will  be  continuous.  We  regard  an  actual  path 
of  the  object  as  measurements  derived  from  a  model  of 
the  predicted  path 


y(^)  =  m(f)-i-X(f) 

where  m{t)  is  a  smooth  function  and  X(t)  is  the  error 
term.  K(f),  m(t)  and  X{t)  are  vector  valued,  i.e., 

X{t)  =  {Xi{t),...,X„{t)f 

m{t)  =  ,  m„(i))^ 

Note  that  the  estimation  of  the  curve  m{t)  is  modeled 
as  a  vector- valued  nonlinear  regression  problem. 

We  desire  to  estimate  the  mean  path  by  m{t).  We 
wish  fh{t)  to  be  “close  to”  the  mean  path  m(t)  and 
“smooth.”  By  “close  to,”  we  mean  that  a  measurement 
of  the  difference  between  rh{t)  and  the  mean  path  (the 
“lack  of  fit”)  is  small.  By  “smooth,”  we  mean  that 
each  of  the  m,(t)  is  “smooth”  in  the  sense  of  belong¬ 
ing  to  a  Sobolov  space.  Thus,  we  require  m(f)  to  be 
a  component-wise  Sobolov  function.  Because  “close¬ 
ness”  and  “smoothness”  are  competing  qualities,  we  se¬ 
lect  m(f)  to  be  that  component- wise  Sobolov  function 
which  minimizes  a  weighted  sum  of  a  lack-of-fit  term 
and  a  smoothing  term.  This  is  the  general  formulation 
for  the  classic  penalized  smoothing  spline  (see  Wegman 
and  Wright,  1983).  The  difference  between  our  formula¬ 
tion  here  and  that  discussed  in  Wegman  and  Wright  is 


that  we  have  a  vector-valued  spline  and  in  the  present 
problem,  we  have  in  mind  a  more  complex  error  struc¬ 
ture  than  simple  i.i.d.  random  variables. 

We  assume  that  the  error  term  behaves  in  accordance 
with  Brownian  motion:  that  is,  1)  the  error  changes  con¬ 
tinuously  over  time  and  the  error  at  any  particular  in¬ 
stant  of  time  appears  to  be  normally  distributed,  2)  the 
error  over  any  interval  depends  on  the  length  of  the  in¬ 
terval  rather  than  on  the  placement  of  the  interval  itself; 
and  3)  the  error  in  any  two  disjoint  intervals  is  inde¬ 
pendent.  We  also  assume:  4)  that  the  expected  error 
at  any  time  t  is  zero.  Finally,  we  assume:  5)  that  the 
error  term  has  greater  variability  between  each  adjacent 
pair  of  joints  than  at  the  joints  themselves.  The  latter 
assumption  is  crucial.  For  example,  an  aircraft  typically 
flies  between  waypoints  on  its  route  of  flight.  A  goal 
of  the  pilot,  perhaps  aided  by  cockpit  automation,  is 
to  cross  the  waypoints.  However,  between  two  consec¬ 
utive  waypoints,  the  aircraft  can  vary  somewhat  from 
a  straight  line,  due  to  errors  in  ground-  or  space-based 
navigation  aids,  on-board  navigation  equipment,  and  pi¬ 
lotage.  Thus,  while  the  object  may  “wander”  somewhat 
from  the  planned  path  on  the  interior  of  each  line  seg¬ 
ment,  the  object  is  expected  to  return  “close  to”  the 
planned  path  at  each  joint. 

More  specifically,  we  assume  that  the  error  term  be¬ 
haves  as  if  it  were  a  Brownian  bridge  between  the  joints 
ti  £uid  ti^i,  that  is,  that  X(t)  is  a  Brownian  motion 
process  conditioned  so  that  ^(^»)  =  =  0.  The 

covariance  structure  of  X  between  the  joints  ti  and 
is  of  the  form 

cov[X]  =  B  0  S 


where  0  is  the  left  Kronecker  product  and  B  is  a  positive 
definite  matrix.  We  wish  to  find  the  covariance  structure 
of  X  between  its  first  and  last  joints  and  tm-  Because 
the  error  term  is  a  Brownian  motion  process,  its  behav¬ 
ior  between  two  adjacent  joints  is  independent  of  its  be¬ 
havior  between  any  two  other  adjacent  joints  (i.e.,  the 
process  has  orthogonal  increments).  Thus  the  covariance 
structure 


of  X  between  the  first  and  last  joints  ti  and  tm  is  the 
block  diagonal  matrix 


/  cov[Xi,i]  0 


0 


cov[Xi^i] 


0 


\ 


0 

0  COv[X,r,,m]  / 


where  each  cov[Xi^i]  is  defined  by 


V  cot.[.Y(<r‘),^(*})]  •••  cov[Xitf<),Xitr‘)] 


for  <  f J  <  . . .  <  <  fj+i  .  Now  cov\Xi^j]  =  0 

when  *  ^  i  because  of  the  independent  behavior.  FVom 
the  standard  theory  of  conditioned  Brownian  motion,  we 
know  that  each  cov\Xi^^  is  of  the  form 


cov[xi^i] 


/  2  fi)(^»+l  ~  fj)  \ 

^  *■  ti+i  -  ti  > 


for  ti  <tf  <t\  <  t,+i.  Letting 


ti+i  -  ti 


B  =  diag{Bi^i) 


and 


t  =  diag{ti^i) 

with  I  being  the  in  x  in  identity  matrix,  we  can  write 
cov[X]  =  B  0  E 


as  the  covariance  matrix  cov[X]  between  the  first  and 
last  joints  ti  and  tm  where  B  is  a  positive  definite  matrix. 
Taking  advantage  of  the  covariance  structure  between 
the  joints,  a  natural  lack-of-fit  term  is 

^  -  m(ti)fi:-\Yj  -  mitj)), 

t=l  j=l 

where  Y  =  and  T 

denotes  transpose. 

Since  m(^)  is  vector- valued, 


J{L{m{t))f{L{mit)))dt=  j  '^{L{rni{t))fdt 


is  appropriate  as  a  smoothing  term. 

Thus,  we  formulate  the  following  minimization  prob¬ 
lem:  Choose  m{t)  to  minimize 


<=i  i=i 


m{tj)) 


0 


4-  u  J  {L{m{t)))'^  {L{m{t)))dt 

where  B  and  E  are  as  described  and  u  >  0,  Miller  and 
Wegman  (1987)  prove  that  this  minimization  problem 
can  be  solved  componentwise,  producing  a  vector  spline 
which  is  a  vector  of  scalar  splines.  They  show  that  the 
objective  function  may  be  rewritten  as 

E  -  z,{tk)) 

f=i  j=i  k-i 

+  pj  {L{Zi{t))fdt) 

where  E  =  PAP^  is  the  spectral  decomposition  (i.e.,  P 
is  orthogonal  and  A  =  diag{Xi)  is  a  diagonal  matrix), 
Wife  =  P^Yk^  and  Z{t)  =  P^m(t).  Since  each  of  the 
p  components  is  non-negative,  the  sum  is  minimized  by 
minimizing  each  of  the  components.  Note  that  there  are 
n  observations  in  p-space. 

3.  Confidence  Bands 

The  confidence  region  is  an  envelope  around  the  pre¬ 
dicted  path.  The  actual  path  of  the  object  is  expected 
to  be  located  within  the  confidence  region  with  a  high 
degree  of  probability.  We  are  given  n  observations  fitting 
the  model  Y (t)  =  m{t)  -f  X(f).  Because  a  vector- valued 
minimization  problem  can  be  solved  component- wise,  we 
can  assume  here  that  T(t),  m(t),  and  X{t)  are  scalar 
valued.  We  assume  that  the  error  term  Jf(f)  is  a  con¬ 
ditioned  Brownian  motion  process  (a  Brownian  bridge) 
between  each  pair  of  points  between  [a,  6];  that  is,  for 
each  t  between  a  and  6,  X{t)  is  normal  with  mean 

=  A  +  (B-yl)L^ 

for  some  real  A  and  B,  variance 

2  2  (*»-<)(< -a) 

(6-a) 

for  some  <7*^,  and  covariance 

cot;[X(s),  X(0]  = 

for  a  <  s  <t  <b.  In  particular,  the  covariance  structure 
of  X{t)  is  of  the  form 

cov[X{t)]  =  B  <g)  E 


where  (g)  is  the  left  Kronecker  product,  B.is  a  positive 
definite  matrix,  and  B"*^  = 

The  solution  to  the  minimization  problem  minimize: 

f=l  j=l  ‘'a 

which  takes  advantage  of  the  coveiriance  structure  of 
X{t)  is  a  cubic  spline  s(x).  We  propose  that  a  (1  —  a) 
100  percent  confidence  region  for  /  is 

C{a,x)  =  s{x)±Za/2y/^ 

where  is  the  a/2  point  of  the  standard  normal  dis¬ 
tribution  and  Vx  is  a  data-based  estimate  of  the  error  at 

X, 

4.  An  Air  Traffic  Control 
Example 

We  apply  the  theory  to  the  modeling  of  the  path  of 
an  aircraft.  While  the  pilot  of  an  aircraft  might  in¬ 
tend  to  fly  directly  from  one  point  to  the  next,  there 
are  many  sources  of  error  preventing  this  direct  flight. 
When  the  pilot  is  using  an  area  navigation  (RNAV)  sys¬ 
tem  designed  for  use  in  the  U.  S.  National  Airspace  Sys¬ 
tem,  the  sources  of  navigation  system  horizontal  error 
are  assumed  to  include  ground  VOR  radiated  signal,  air¬ 
borne  VOR  receiver  equipment,  area  navigation  equip¬ 
ment,  and  pilotage.  (See  DOT/FAA,  1975.)  In  addition, 
it  is  assumed  that  the  errors  attributable  to  these  four 
sources  are  independent  normal  distributions  that  may 
be  combined  root-sum-square  (RSS)  fashion.  (See  also 
ICAO,  1985.)  The  2s  (95th  percentile)  values  of  the  er¬ 
ror  components  assumed  in  the  U.S.  are  summarized  in 
Table  1. 

Table  1.  Horizontal  Error  Components  of  the 
VOR/DME  System 


Error  Component 

2s  Error 

VOR  Equipment 

Ground  VOR  radiated  signal 

1.2° 

Airborne  VOR  receiver 

2.7° 

Course  setting  equipment 

1.6° 

Pilotage 

2.3° 

Suppose  that  the  aircraft  is  navigating  by  reference 
to  the  ground-based  navigation  aids  (NAVAIDs)  of  the 
air  traffic  control  (ATC)  system  (see  FAA,  1992).  When 
the  route  between  two  NAVAIDs  is  not  part  of  an  es¬ 
tablished  airway  or  route  and  is  below  flight  level  600 


(approximately  60,000  feet  in  altitude),  then  air  traf¬ 
fic  controllers  must  protect  the  airspace  along  the  route 
for  four  miles  on  each  side  of  the  route  imtil  51  miles 
from  the  NAVAID  and  then  increasing  along  a  4.5  de¬ 
gree  angle  to  ten  miles  on  each  side  at  130  miles  from  the 
NAVAID.  See  Figure  1.  Thus,  we  can  begin  to  estimate 
the  allowable  displacement  along  the  route. 

The  shape  of  the  protected  route  width  suggests  that 
the  allowable  displacement  could  be  defined  as  a  Brow¬ 
nian  motion  process  conditioned  to  be  (essentially)  zero 
at  both  ends  of  the  route  segment,  i.e.,  at  the  NAVAIDs. 
If  the  displacement  at  time  t  is  represented  by  V(f), 
then  Y  (f )  is  a  Browniein  motion  process  conditioned  by 
y(fi)  =  0  =  ^(^2),  when  Y{ti)  represents  the  displace¬ 
ment  at  time  #1  at  the  first  NAVAID  and  Y (*2 )  represents 
the  displacement  at  time  f2  at  the  second  NAVAID. 


< - 2e0nrni  - > 


Figure  1.  Protected  Route  Width  (ATC) 

We  know  that 

Y{t)  ~ 

h  “  ti 

Letting  =  0  and  t2  =  260,  the  v^lriance  is  of  the 
form  (7^(260  — f)f/260,  showing  that  the  largest  displace¬ 
ment  is  possible  at  the  midpoint  {t  =  130)  between 
the  two  NAVAIDs.  We  estimate  a  by  assuming  that 
1.96^^0-2(260  —  f)f/260  represents  the  95th  percentile  of 
the  distribution  and  solving  the  equation 

10  =  1.96  (260  -  130)  ■  130/260 

for  a  .  The  solution  is  tr  ~  0.632830279.  Figure  2  com¬ 
pares  two  displacement  regions  on  one  side  of  the  route: 
the  ATC-defined  region  and  the  95th  percentile  of  the 


derived  distribution. 


Figure  2.  Protected  Route  Width  (Conditioned  Brow¬ 
nian  Motion  Process) 


5.  Conclusions 

In  this  paper  we  outline  a  strategy  for  the  estimation 
of  the  path  of  a  vehicle  and  a  confidence  region  around 
that  path.  For  an  important  class  of  vehicles  which  in¬ 
cludes  commercial  airliners,  the  path  is  a  curve  in  three- 
dimensional  space.  This  places  path  modeling  in  the 
domain  of  differential  geometry  because  we  wish  a  path 
to  be  a  curve  in  three-space  with  at  least  two  continuous 
derivatives  (namely  velocity  and  acceleration).  A  path  is 
modeled  as  a  vector- valued  function  with  two  continuous 
derivatives.  We  can  assume  therefore  that  the  path  is  a 
curve,  an  element  of  a  vector- valued  Sobolov  space.  An 
observed  path  is  modeled  as  the  curve  with  added  noise. 

The  nature  of  the  additive  noise  departs  from  the  typ¬ 
ical  white  Gaussian  noise  assumption  both  because  of  a 
likely  correlation  in  the  noise  structure  and  also  because 
the  nature  of  motion  towards  a  destination  implies  less 
veiriability  at  the  origin  and  destination  of  the  vehicle 
motion.  For  this  reason  we  assume  a  Brownian  bridge 
model  for  the  additive  noise.  The  path  being  modeled 
as  a  vector- valued  function  with  two  continuous  deriva¬ 
tives  taken  together  with  a  Brownian  bridge  error  struc¬ 
ture  implied  a  fit  to  the  path  that  is  a  vector  of  cubic 
splines  albeit  not  a  traditional  cubic  spline.  Miller  and 
Wegman  (1987)  develop  a  theory  of  vector  splines  and 
develop  computational  algorithms  for  certain  classes  of 
vector  splines.  These  methods  apply  to  the  setting  in 
which  the  error  structure  is  a  Brownian  bridge. 


We  conclude  the  paper  with  a  comparison  of  the  pro¬ 
tected  route  width  of  the  U.S.  air  traffic  control  (ATC) 
system  and  the  estimated  95%  confidence  band  assuming 
the  spline-Brownian  bridge  model.  If  the  model  is  valid 
for  aircraft  movement,  then  it  can  be  suggested  that 
current  air  traffic  control  procedures  reserve  too  much 
airspace  near  navigation  aids  and  too  little  airspace  mid¬ 
route.  Thus  the  qualitative  shape  difference  is  of  con¬ 
cern,  and  our  work  suggests  that  a  different  model  for 
protected  airspace  may  be  appropriate.  We  note  that 
we  predicated  our  confidence  band  on  matching  maxi¬ 
mum  route  widths  with  the  ATC  protected  route  width 
specification.  Because  data  to  actually  calibrate  a  95% 
confidence  band  is  not  easily  available,  conclusions  about 
the  adequacy  of  the  protected  route  width  in  quantita¬ 
tive  terms  are  not  warranted.  A  realistic  empirical  study 
of  actual  aircraft  movement  along  ATC  routes  would 
be  fruitful  since  it  is  possible  that  an  accurately  cali¬ 
brated  estimate  of  an  aircraft  path  and  its  confidence 
region  could  conceivably  remove  excessively  restrictive 
protected  route  width  specifications.  In  addition,  the 
confidence  region  parameter  values  used  in  the  example 
are  derived  from  nominal  values  for  navigation  accuracy 
used  for  all  navigation  systems  as  a  maximum  error  bud¬ 
get.  However,  individual  navigation  systems  are  capable 
of  less  error.  The  nominal  values,  being  conservative, 
require  that  more  airspace  be  reserved  for  each  aircraft. 
Using  error  values  tailored  to  individual  navigation  sys¬ 
tems  might  also  allow  better  utilization  of  airspace.  Far¬ 
ther  studies  could  determine  (a  library  of)  error  values 
for  the  confidence  region  parameters. 
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